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Abstract 


To understand further the thermal abuse behavior of large format Li-ion batteries for automotive applications, the one-dimensional modeling 
approach formulated by Hatchard et al. [T.D. Hatchard, D.D. MacNeil, A. Basu, J.R. Dahn, J. Electrochem. Soc. 148(7) (2001) A755—A761] was 
reproduced. Then it was extended to three dimensions so we could consider the geometrical features, which are critical in large cells for automotive 
applications. The three-dimensional model captures the shapes and dimensions of cell components and the spatial distributions of materials and 
temperatures, and is used to simulate oven tests, and to determine how a local hot spot can propagate through the cell. In simulations of oven abuse 
testing of cells with cobalt oxide cathode and graphite anode with standard LiPF, electrolyte, the three-dimensional model predicts that thermal 
runaway will occur sooner or later than the lumped model, depending on the size of the cell. The model results showed that smaller cells reject 
heat faster than larger cells; this may prevent them from going into thermal runaway under identical abuse conditions. In simulations of local hot 
spots inside a large cylindrical cell, the three-dimensional model predicts that the reactions initially propagate in the azimuthal and longitudinal 


directions to form a hollow cylinder-shaped reaction zone. 
© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


Small lithium-ion (Li-ion) cells and multistring battery packs 
are now commonly used in consumer products. Larger packs 
(>200 Wh) are limited to niche applications like satellites, in 
large part because of safety concerns. With increasing interest 
in large-format Li-ion batteries for automotive applications, we 
need to better understand the abuse tolerance of these batteries. 

Li-ion batteries present an attractive energy storage device for 
electric vehicles (EVs) and hybrid electric vehicles (HEVs) for 
several reasons. The high cost of nickel may give Li-ion batteries 
a cost advantage over nickel/metal-hydride (Ni/MH) batteries. 
Ni/MH cells, which are widely used in HEVs, have a voltage 
of about 1.2 V, whereas Li-ion cells typically have a voltage of 
about 3.6 V. So the same system voltage can be achieved with 
one-third the number of Li-ion cells as Ni/MH cells. Reducing 
the number of cells dramatically increases reliability because 
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the cells are usually connected in series. Thus, the failure of a 
single cell could deactivate an entire module. Li-ion cells have 
up to twice the gravimetric energy density and two to three times 
the gravimetric power density of Ni/MH cells, which improve 
the energy efficiency of HEVs. Plug-in hybrid electric vehicle 
(PHEV) or EV applications, which require higher energy/power 
battery systems, may benefit even more. Lastly, recent evidence 
suggests that Li-ion batteries have the potential of achieving bet- 
ter low-temperature performance than Ni/MH batteries. These 
advantages make Li-ion an attractive candidate for advanced 
automotive energy storage systems, but adoption of the tech- 
nology is delayed by problems of cost, life, low-temperature 
performance, and abuse tolerance. 

Abuse tolerance is an especially daunting problem as it relates 
to safety issues. Acceptable abuse tolerance in small consumer 
Li-ion cells (<3 Ah) and packs (<150 Wh) has been achieved by 
the generous (and often redundant) use of safety devices (such 
as shutdown separator, positive temperature coefficient [PTC], 
current interrupt, electronic, and pressure vent devices) [1,2]. 

However, achieving good abuse tolerance with the larger Li- 
ion cells (>6 Ah) and packs (>200 Wh) needed for automotive 
applications remains challenging. The high-power requirements 
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Fig. 1. Comparison of chemistry box model for three-dimensional model with results of Hatchard et al. [5] for oven heating (155 °C) of LiCoO2/graphite 18650 cell. 
Figure at left is from the box model for three-dimensional model developed here and figure at right is taken from Hatchard. 


of the HEV/PHEV/EV applications prevent the use of PTC and 
other safety devices. Larger cells are intrinsically more vulnera- 
ble to thermal runaway because they have higher energy content. 
Increasing cell size lowers the surface:volume ratio, which in 
turn reduces cooling area per volumetric heat generation. In 


addition, increasing the size of a cell increases the temperature 
difference inside a cell for two reasons: higher heat flux caused 
by smaller surface:volume ratio, and higher thermal resistance 
caused by a larger heat diffusion distance. A wider cell tem- 
perature distribution for a given surface heat transfer condition 


Table 1 

Physical and kinetic parameters used for abuse simulations 

Symbol Description Value Refs. 
Asgi SEI-decomposition frequency factor 1.667E15 (s7!) [5] 
Ane Negative-solvent frequency factor 2.5E13 (s7!) [5] 
Ape Positive-solvent frequency factor 6.667E13 (s7!) [5] 
Ae Electrolyte decomposition frequency factor 5.14E25 (s7!) [6] 
Ea sei SEI-decomposition activation energy 1.3508E5 (J mol!) [5] 
Fane Negative-solvent activation energy 13508E5 (J mol!) [5] 
Eape Positive-solvent activation energy 1396E5 (Jmol!) 

Eze Electrolyte decomposition activation energy 2.74E5 (J mol!) [6] 
Csei0 Initial value of Csei 0.15 [5] 
Cneg0 Initial value of Cneg 0.75 [5] 
ao Initial value of œ 0.04 [5] 
CeO Initial value of ce 1 [5] 
Mgei Reaction order for Csei (Eq. (7)) 1 [5] 
Mne,n Reaction order for Cneg (Eqs. (10)-(11)) 1 [5] 
Mpe,p! Reaction order for «œ (Eq. (14)) 1 [5] 
Mpe,p2 Reaction order for (1 — a) (Eq. (14)) 1 [5] 
Me Reaction order for ce (Eq. (17)) 1 

tseid Initial value of tsei 0.033 [5] 
Hei SEI-decomposition heat release 257 (J g7!) [5,6] 
Ane Negative-solvent heat release 1714 (J g7!) [5,6] 
Ape Positive-solvent heat release 314 (J gl) [5,6] 
Hele Electrolyte decomposition heat release 155 (J g7!) [6,20] 
We Specific carbon content in jellyroll 6.104E5 (g m~’) 

Wp Specific positive active content in jellyroll 1.221E6 (g m7’) 

We Specific electrolyte content in jellyroll 4.069E5 (g m~?) 


Volumetric heat capacity of jellyroll 


2.789E6 (J m~? K7!) 


478 G.-H. Kim et al. / Journal of Power Sources 170 (2007) 476-489 


raises the possibility of local onset temperatures that trigger cell 
thermal runaway. 

A more comprehensive understanding about the thermal 
behaviors of Li-ion batteries at elevated temperatures is required 
to design a large, abuse-tolerant cell. The U.S. Department of 
Energy is supporting research and development at national lab- 
oratories (ATD and BATT Programs) and FreedomCAR battery 
developers (U.S. Advanced Battery Consortium Programs) to 
investigate methods and materials that address the abuse tol- 
erance of Li-ion cells and packs. The development of small, 
abuse-tolerant Li-ion cells for consumer use relies heavily on 
“build and break” cycles. For example, shutdown separator spec- 
ifications were developed and verified based on abuse testing 
large numbers of cells. However, for large automotive cells, this 
trial-and-error design process does not seem viable, because 
it is expensive and time-consuming for large cells; moreover, 
the build and break test would provide only limited infor- 
mation for effectively improving the abuse tolerance of the 
cells. 

The objectives of this study are to provide better insight about 
the behaviors of cells under abuse conditions by developing a 
three-dimensional thermal abuse model for Li-ion batteries, and 
to accelerate the development of abuse-tolerant Li-ion batteries 
for vehicle applications by establishing a tool and methodology 
to support the design. 


2. Model development 


Designing large Li-ion batteries has raised questions such as: 
How do the cell sizes and shapes affect abuse tolerance? How 
important is heat transfer from the tabs to good abuse tolerance? 
Are shutdown separators effective in large cells? If so, what 
is the optimal temperature for onset of shutdown? Will thermal 
runaway propagate to other cells? Under which scenarios? What 
are the safety margins for key design parameters? 

A multidisciplinary science and engineering approach is 
required to understand the issues of Li-ion battery abuse behav- 
iors, which are related to the various aspects from material 
properties to assembly methods. However, the key to resolving 
the abuse tolerance problem in Li-ion batteries lies in thermal 
engineering—controlling the heat generation and the internal 
and external heat transfer. 

Researchers have tried to model the abuse reactions that 
occur in Li-ion batteries [3-8]. Hallaj et al. [8] developed a 
one-dimensional thermal model and identified conditions under 
which large cells could go into thermal runaway. Botte et al. 
[3] modified a DUAL model [9] by adding a term to the energy 
balance to account for an exothermic chemical reaction at the 
negative; they found that the chemical reaction of the nega- 
tive was negligible under normal conditions. Richard and Dahn 
[10] showed how accelerating rate calorimetry (ARC) studies 
of small samples could be used to predict oven abuse results 
for full cells. Hatchard et al. [4] showed the importance of con- 
sidering radiative heat transfer in oven abuse testing. Hatchard 
et al. [5] later presented a model for predicting oven abuse 
tests based on expressions for the exotherms of the electrodes 
obtained by ARC and differential scanning calorimetry studies. 


This approach was extended by Spotnitz and Franklin [6], who 
considered an array of possible exothermic reactions and pre- 
sented approaches for simulating a variety of abuse tests (oven, 
short-circuit, overcharge, nail, crush); this paper also contains 
an extensive review of experimental results. Yamauchi et al. 
[7] described an approach for simulating short-circuit behavior. 
Recently, Spotnitz et al. [11] described an approach for esti- 
mating the stability of packs when one cell goes into thermal 
runaway. 

Abuse thermal behaviors of Li-ion batteries, especially large 
ones, have intrinsically three-dimensional features. Chen et al. 
presented three-dimensional thermal models for prismatic [12] 
and cylindrical [13] cells; these papers show that asymmetric 
temperature distributions arise from geometric effects and from 
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Fig. 2. Simulation of effect of oven temperature on thermal response of 18650 
cell: (a) temperature vs. time; (b) heating rate terms for 150 °C oven; (c) heating 
rate terms for 155°C oven. 
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convective heat transfer. Pesaran et al. [14] and Bharathan et al. 
[15] demonstrated the utility of three-dimensional models for 
improving the current and temperature distributions in batteries. 
However, these papers [12-15] did not consider abuse behaviors. 
Thermal abuse behaviors of Li-ion batteries are greatly affected 
by the local distributions of heat and materials. Therefore, 
three-dimensional modeling for capturing geometrical thermal 
effects on battery thermal abuse behavior is required. A three- 
dimensional model can capture the thermal paths inside a cell 
by addressing the geometries and properties of cell components 
[16]. Various thermal boundary conditions can be simulated 
for considering local cooling or heating effects. In addition, a 
three-dimensional analysis can address the effects of nonuni- 
form distributions on heat release from the component reactions 
in a high-temperature environment. 

To develop a three-dimensional thermal abuse model, we 
listed possible exothermic component reactions, and then imple- 
mented the abuse chemistry into the three-dimensional battery 
model. In this study, a commercial, finite-volume method 
software known as FLUENT® was used to develop the three- 
dimensional thermal abuse model for Li-ion batteries. The 
thermal abuse model does not explicitly consider internal short 
circuits that might occur during oven testing due to soften- 
ing of the separator. However, localized heating of the cell is 
considered. 


2.1. Energy balance 


Thermal conduction is assumed to dominate internal heat 
transport in the system, which is true while the convective heat 
transfer through the net flow of electrolyte and production gas 
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is limited. The temporal change of stored thermal energy in 
arbitrary control volume in the system equals to the difference 
between heat generation and net heat diffusion out, as shown in 
the following equation: 


ð T 
Pe L LYEVT) +S, 
S = Sabuse_chem + Sjoul + Scombustion + ` + (1) 


where p (g cm~?) is the density, Cp J g7! K7!) the heat capacity, 
T (K) the temperature, t (s) the time, k (W cm~! K7!) the thermal 
conductivity, and S (W cm~?) represents the heat source/sink 
terms. The volumetric heat generation that originated from Li- 
ion battery component reactions at increasing temperatures was 
listed in Eq. (2) and included in the three-dimensional battery 
thermal analysis: 


Sabuse_chem = Ssei + Sne + Spe + Sele + Snb (2) 


where Sse is the heat from the SEI decomposition reaction, She 
the reaction between the negative active material and electrolyte, 
Spe the reaction between the positive active material and elec- 
trolyte, Sele the electrolyte decomposition, and Spp is the reaction 
between the negative active and binder. 

One difficulty for the three-dimensional modeling of batteries 
is that extensive computational resources are needed to resolve 
the small-scale structures such as layered winding of electrodes, 
separators, and current collectors in a reasonable grid size. In 
this study, sub-grid scale structures are modeled as a continuous 
material with non-isotropic properties that depend on the layered 
directions of sub-grid structures. Orthotropic properties such as 
thermal conductivity of layered structure can be measured val- 
ues, or can be evaluated with flux conservation relations through 


=> 
o 
~ 

o o o 

> D œ = 


Degree of Reaction 


o 
N 


oO 


0 30 60 90 120 150 180 


Time (min) 


1 


~ 
Q 
= 


o 
œ 


o 
D 


o 
a 


Degree of Reaction 


o 
N 


0 30 60 90 120 150 180 
Time (min) 


Fig. 3. Computed heating rate from the component reactions in simulated oven tests. 150°C: (a) component terms, (b) degree of reaction; 155 °C: (c) component 


terms, (d) degree of reaction. 
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the layers and along the layers of thickness l; (Eqs. (3) and 
(4): 


e Layer-normal direction: 


ER Dili 
kn = > i/ki (3) 


e Layer-parallel direction: 


‘Kil; 
ee ao (4) 


The heat exchange between the system and the ambient is 
computed through the thermal boundary conditions assigned 
at the computational domain boundaries. Combined thermal 
boundary conditions are applied to consider both the radiation 
and the convection heat transfer contributions at the cell surface. 
Convective heat flux out to ambient is evaluated in the following 
equation: 


oe = h(Tsurt — Tamb) (5) 


where h is a convection heat transfer coefficient. For oven test 
simulations, Eq. (6) is used for radiative heat flux out to ambient 
from cell surface: 


4 4 
radi = €0(Tyur¢ — Tamb) (6) 


where €e is the emissivity of the cell surface and ø is the 
Stefan—Boltzmann constant. The oven chamber is assumed to 
be an isothermal cavity where the irradiation corresponds to 
emission from a blackbody, and the cell surfaces are considered 
as gray diffuse surfaces. 


2.2. Thermal abuse reaction models 


Li-ion battery thermal abuse reactions are modeled for three- 
dimensional battery simulation by listing component reactions 
that occur at elevated temperatures. Reactions involving lithium 
metal (which would be important in an overcharge test [6]) and 
combustion reactions are not considered in the current study, 
but the approach presented here could readily accommodate 
such reactions. Therefore, our model assumes that the flames 
are successfully prevented for the simulated cells that are not 
overcharged. 


2.2.1. Solid electrolyte interface decomposition reaction 

The negative electrode is protected from direct reaction 
with solvent by an ionically conducting film called the solid 
electrolyte interface (SEI). The layer is meta-stable and can 
decompose exothermically at 90-120 °C [17,18]. Following the 
formulation of Hatchard et al. [5], this reaction can be expressed 
in the following equations: 


Ea sei f 
Rsei(T, Csei) = Asei exp |- a coset (7) 


Ssei = Hsei We Rsei (8) 


= — Rsei (9) 


where Cgej is the dimensionless amount of lithium-containing 
meta-stable species in the SEI. Rgej (s7}), Asei (s71), Ea sei 
J mol`!) are reaction parameters. Specific heat release Hsei 
(J g7!) is evaluated as Joule per gram of carbon, and We (g m~’) 
is the volume-specific carbon content in the jellyroll. 


2.2.2. Negative-solvent reaction 
At elevated temperatures (>120 °C), an exothermic reaction 
between intercalated lithium and electrolyte can occur (Eqs. 
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Fig. 4. Simulated effect of cell size in 150 °C oven test: (a) temperature vs. time; 
(b) component heat terms for D18H65 cell; (c) component terms for D50H90 
cell. 
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Fig. 5. Simulation of effect of cell size and surface:volume ratio (A/Vjr) on different temperature oven storing. (See Table 2, Cell A (A/Vj,=48.57), Cell B 
(A/Vjr =25.01), Cell C (A/Vjr =29.14).): (a) 160°C oven test; (b) 150°C oven test; (c) 140°C oven test. 


(10)-(13)). Again, following the formulation of Hatchard et al. 
[5]: 


tsei Fane 
Rne(T, Ce, Cneg, tsei) =Ane EXP |- Ce exp | — ; 
sei,ref RT 


tunneling regime (10) 
tsei,ref Ea,ne 
Rne(T, Ce, Cneg, tsei) = Ane ( a ) Cae exp j-e] g 
diffusion regime (11) 
She = Ane We Roe (12) 
dcneg 
—=-=-R 13 
J ne (13) 


where Cneg is the dimensionless amount of lithium intercalated 
within the carbon; tsei a dimensionless measure of SEI layer 
thickness that reflects the amount of lithium in the SEI; and 
Rite “Acie >), Eane (J mol— ') are the reaction parameters. 
Specific heat release Hne (J g7!) is evaluated as Joule per gram 


Table 2 
Cell sizes and surface:volume ratios (A/Vjr) used for oven test simulations 

Cell A Cell B Cell C 
Vir (cm?) 10.52 157.1 157.1 
A/V; (m7!) 48.37 25.01 29.14 


of carbon, and W, (g m7?) is the volume-specific carbon content 
in the jellyroll. 


2.2.3. Positive-solvent reaction 

In the oxidized state, the positive material reacts directly with 
the electrolyte. Or, the positive active material can decompose 
exothermically and emit oxygen that can react exothermically 
with the electrolyte. In any case, the chemical reduction of the 
positive active material with the electrolyte is highly exothermic 
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Fig. 6. Schematic of model geometry and a computational mesh. 
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Fig. 7. Three-dimensional simulation of 155°C oven test for a 18650 cell: (a) temperature distribution at 55 min; (b) heat generation for positive-solvent and 


negative-solvent reactions at 55 min. 


(Eqs. (14)-(16)). Following the formulation of Hatchard et al. 
[5]: 


E 
Rpe(T, 0, Ce) = Apea” (1 — oP"? exp |- se| (14) 


RT 
Spe = Hpe Wp Rpe (15) 
da 
dt = Rpe (16) 


where æ is the degree of conversion. Rpe (s7!), Ape (s—}), Ea Pe 
J mol!) are reaction parameters. Specific heat release Hpe 
(J g7!) is evaluated as Joule per gram of positive-active content, 
and Wp (g m`?) is the volume-specific, positive-active content 
in the jellyroll. 


2.2.4. Electrolyte decomposition reaction 
The electrolyte can decompose exothermically at elevated 
temperatures (>200 °C), as expressed in the following equations: 


E 
R(T, Ce) = Ac exp |- sc ce (17) 


R 


Sele = He We Re (18) 
dce 

e —_R 19 
di e (19) 


where ce is the dimensionless concentration of electrolyte and 
Re (s7!), Ae (s7!), Eae (J mol`!) are the reaction parameters. 
Note that ke is a frequency factor rather than a rate constant. 
Specific heat release He (J g7!) is evaluated as Joule per gram 
of electrolyte, and We (g m~?) is the volume-specific electrolyte 
content in the jellyroll. 


2.3. Model verification 


We verified the chemistry portion of the three-dimensional 
implementation by comparing a simple lumped! thermal model 
with the one-dimensional results presented by Hatchard et al. 
[5]. The chemistry portion of the model reproduced the results of 


! The term “lumped model” indicates that internal temperature gradients are 
relatively small so that a single temperature can be used for the entire model 
domain. 
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Hatchard et al. [5], with small differences caused by the lumped 
versus one-dimensional treatment (Fig. 1). 


3. Results 


The model parameters used in this study were collected 
from the literature. Open information is limited, but Table 1 
presents the values of model parameters used to describe the 
LiCoO>/graphite chemistry at full charge. The lumped model 
results were used to investigate the effects of cell size and oven 
temperature on the thermal abuse behavior of a cell under the 
oven test. Three-dimensional oven test simulation results are 
presented to address the impact of temperature nonuniformity 
and material distributions through a cell. The three-dimensional 
model was used to simulate a localized heat release and the 
propagation of the chemical reactions inside a cell. 


3.1. Lumped model-oven test simulation 


Oven test simulations were carried out with the lumped model 
to investigate the impact of oven temperature and cell size. The 
cells and batteries were assumed to be initially at a normal oper- 
ating temperature (35 °C), and then suddenly placed in an oven 
that was preheated to the desired test temperature. The oven tem- 
perature was kept constant during the test. Natural convection 
and radiative heat transfer were considered at the cell surface 
(h=7.17Wm~ K, ¢=0.8). 

We examined the effect of oven temperature on the thermal 
behavior of an 18650-size cell (cylinder, 18 mm diameter, 65 mm 
height) for two oven temperatures; 150 and 155°C. We intro- 
duced thermal mass-specific heating rates, Q/Mcp (°C min” !) 
to quantify the internal heat generation and heat exchange with 
ambient. Fig. 2(a) shows the cell temperature variations for both 
oven temperature cases. Internal heat generation that resulted 
from thermal abuse reactions and heat transferred from ambient 
were plotted with the sum of both contributions for each oven 
temperature in Fig. 2(b) and (c), respectively. The cell placed in 
the 150°C oven did not go into thermal runaway; the cell heating 
in the oven at 155°C did go into thermal runaway. The initial 
heating rate profiles are similar in both oven temperature cases. 
The initial heating rates (t <20 min) were determined mainly by 
the heat transfer from ambient. 

In the 150°C oven heating case, heat generation and heat 
rejection came to balance at about 40 min. The net heating rate 
increased gradually afterward as the reactions slowly went on 
until about 110 min when the reactants were depleted. On the 
other hand, in the 155°C case, the heat generation rate rapidly 
increased to exceed the heat rejection to ambient at around 
50 min and led to thermal runaway. 

To better understand the effect of oven temperature, the heat 
release by each component was closely examined. Heat gener- 
ation from each component reaction was quantified and plotted 
along with the net heat generation for the 150°C oven tem- 
perature case in Fig. 3(a). Time variations of reaction degree 
of the component reactions are shown in Fig. 3(b). Fig. 3(c) 
and (d) presents the same quantities for the 155°C oven tem- 
perature case. In the initial heating period (t< 15 min), the SEI 


decomposition reaction contributed the most to total heat gener- 
ation in both temperature cases to show the first local maximum. 
This heat was released between 80 and 140°C and was mostly 
used to heat the cell. At 20 min <¢< 30 min, the negative-solvent 
reaction led to total heat generation. The reaction slowed as 
the diffusion thickness for reacting Li-ion in the negative parti- 
cles increased. The positive-solvent reaction was autocatalytic 
[19]; that is, the reaction initially proceeded slowly but accel- 
erated rapidly as reaction products that promoted the reaction 
were formed. Positive-solvent reaction heat release increased 
gradually and dominated the total heating rate at t>60 min in 


(a) 250 
200 
150 


100 


Call Temperature (°C) 


= 
o 
= 


SEI Lumped 
— SEI 3D 

Negative-Solvent Lumped 
—— Negative-Solvent 3D 

Positive-Solvent Lumped 
—— Positive-Solvent 3D 


Specific Heat Generation (°C/min) 


= 
O 
~ 


Nondimensional Amount 


Time (min) 


Fig. 8. Comparison of three-dimensional and lumped model results for 155°C 
oven test of a 18650 cell: (a) temperature; (b) heat generation of component 
reactions; (c) reaction quantity. 
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Fig. 9. Schematic of temperature gradients computed from lumped and three- 
dimensional models. 


the 150°C oven heating case, and at t>40 min in the 155°C 
oven heating case. In the 155°C case, a sudden temperature 
increase (at about 60 min) caused by the positive-solvent reac- 
tion enhanced the negative-solvent reaction again and triggered 
the electrolyte decomposition. But the reaction ceased after cool- 
down to ambient temperature. In the 150°C case, internal heat 
generation occurred slowly enough for the generated heat to be 
rejected to ambient without thermal runaway. 
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The effect of a cell’s size on its thermal behavior was inves- 
tigated for a fixed oven temperature. A small cell (D18H65, 
18 mm diameter and 65 mm height) and a large cell (D50H90, 
50mm diameter and 90 mm height) were simulated for 150°C 
oven temperature heating. Fig. 4(a) shows the cell temperature 
variations for the both size cell cases. Internal heat genera- 
tion caused by thermal abuse reactions and heat transferred 
from ambient was plotted with the sum of both contributions 
for each size of cell, as shown in Fig. 4(b) and (c), respe- 
ctively. 

The rate of temperature rise was initially much slower in a 
large cell than in a small cell because the volume-specific heat 
exchange area was smaller in the large cell. So, the small cell 
heated rapidly but stayed at slightly higher temperature than 
the oven temperature while the exothermic reactions slowly 
occurred. In the small cell case, heat rejection to ambient bal- 
anced with heat generation so thermal runaway was averted. 
In the larger cell (D50H90), heat generation started slowly but 
increased continuously. Heat rejection started much later (about 
50 min) in the larger cell and the magnitude of rejection heat 
was relatively small compared to heat generation, which led to 
thermal runaway. 


negative-solvent reaction 


Fig. 10. Simulation of 155°C oven test for D50H90 cell using three-dimensional model. 
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Heat generation rates per cell thermal mass are plotted in 
Fig. 5 for the oven heating time at three oven temperatures with 
three cells. Table 2 presents the jellyroll volumes and volume- 
specific heat transfer areas for the tested cells. Cell A and Cell B 
represent D18H65 and D50H90 size cylindrical cells, respec- 
tively; Cell C represents a cell that has same volume with 
D50H90 cell but a different form factor to have a larger heat 
exchange surface. The number for the volume-specific surface 
area of Cell C was calculated from an elliptical cylinder cell. 

Heat exchange surface area per jellyroll volume, A/Vj;, is a 
critical factor for determining initial heating and heat dissipa- 
tion to ambient. On the other hand, the heat generation rate is 
sensitive to the cell temperature. For larger A/Vj;, the cell leads 
to thermal runaway faster at high oven temperature (160 °C), but 
the heat rejection is large enough to balance with the heat gener- 
ation at low oven temperature (140°C) so thermal runaway can 
be avoided. The result implies that increasing A/V; lowers the 
onset oven temperature, which triggers the cell thermal runaway, 
and large A/V;, shortens the time for oven heating to make the 
cell undergo thermal runaway, if it occurs. 


3.2. Three-dimensional, lithium-ion thermal abuse model 


We used the model based on finite-volume methods to 
conduct three-dimensional thermal abuse simulations for Li- 
ion batteries. The parameters for atypical fully charged 
LiCoO2/graphite chemistry were used (see Table 1). Fig. 6 
shows a schematic of a sample model geometry and compu- 
tational mesh used in the simulations. One-half model with 
symmetry surface was used. The total number of computa- 
tional cells in the meshes is about 42,000. The computational 
cells are hexahedral and ranged in volume from 3.13E—10 to 
5.44E—9 m°. The model geometry includes the cell compo- 
nent details such as aluminum/copper current collectors and 
terminals, insulator blocks, aluminum container, and cell core 
winding. We addressed the material property variations in the 
continuous calculation domain by dividing the domain into 
several computational zones. For the cells studied here, the 
container can is electrically and thermally connected to one of 
the terminals (Fig. 6). NREL’s Linux cluster was used for the 
three-dimensional computations. 


3.2.1. Three-dimensional oven test simulation 

Even though the oven test is not a conspicuous case in which 
the three-dimensional effects are critical, 155 °C oven test cases 
for the different size cells were simulated three dimension- 
ally, and the results were compared with the lumped analysis 
results to see the impact of nonuniform distributions. The model 
incorporates the heat released from the component reactions, 
SEI decomposition, negative-solvent reaction, positive-solvent 
reaction, and electrolyte decomposition reaction. The thermal 
boundary condition for the cell exterior surface addresses natu- 
ral convection and radiation heat exchange (1 =7.17 Wm?K, 
€=0.8). Black-body irradiation and gray surface conditions 
were assumed. 

Fig. 7 presents the contours of temperature and heat gen- 
eration for the D18H65 small size cell oven test after storing 


at 155°C for 3300s. At the specified time, the cell undergoes 
thermal runaway. The highest temperature appears at the cell 
center where the chemical reactions are stronger. The tempera- 
ture gradient is larger in the layer normal (radial) direction. A 
larger temperature gradient does not directly mean larger heat 
flux because of the nonisotropic thermal conductivity of the cell 
core. The bottom of the cell core winding is thermally connected 
to the can through the current collector. So, the lower section of 
the core is cooled better and results in lower temperature and 
heat generation than the upper section. 
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Fig. 11. Comparison of lumped and three-dimensional models for 155°C oven 
test of a D50H90 cell: (a) temperature; (b) heat generation of component reac- 
tions; (c) reaction quantity. 
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The temporal variations of three-dimensional model cell 
average values for temperature, heat generation, and reaction 
parameters were compared with the lumped analysis quantity 
variations in Fig. 8 for the D18H65 cell. Nonuniform distribu- 
tion affects the results in two aspects: (1) heat exchange between 
the cell and ambient environment is more accurately calculated 
by capturing the local temperature profiles at the surface and 
(2) the spatial variation of temperature and the local reaction 
parameters are captured so the chemical reaction rates that are 
susceptive to local variables can be evaluated more accurately. 
As shown in Fig. 8, the three-dimensional simulation results are 
generally matched with the lumped system analysis, but the start 
of thermal runaway is moved up by several minutes. Because of 
the (relatively) large A/V;, of D18H65, the cell heats rapidly 
and starts to reject heat to ambient at about 20 min. During the 
heat rejection period, the captured temperature profiles in the 
three-dimensional model make the heat rejection rates at the 
cell’s surface be lower, and reaction rates at the cell’s inner core 
region are higher than in the lumped system analysis (Fig. 9). 
So, the thermal runaway is advanced by several minutes in the 
three-dimensional simulation compared to the lumped analysis 
case. 

Fig. 10 presents the contours of temperature and heat gener- 
ation for the D50H90 size cell oven test after storage at 155°C 
for 3860s. At the specified time, the cell undergoes thermal 
runaway. Spatial distribution nonuniformity is more apparent 
than with the D18H65 case. The highest temperature appears 
at core center region. The reactants of positive-solvent reaction 
are already depleted at the center region of the cell so that the 
reaction occurs only in the outer rim and at both ends (near 
the top and bottom) of the winding at the specified time. Since 
the heat generation rate is much larger than the surface heat 
rejection rate, temperature distribution is dominated by the heat 
generation rather than by the location of heat transfer surfaces. 
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We examined the nonuniform distribution effects in the larger 
cell oven test case by comparing the averaged values from the 
three-dimensional model results with the lumped model results 
(Fig. 11). The three-dimensional simulation results generally 
matched with the lumped system analysis. But, the start of ther- 
mal runaway was delayed by a few minutes and the maximum 
temperature was higher. Due to the (relatively) small A/V; of 
D50H90, the cell heated slowly and did not reject heat to ambient 
until about 50 min. (Thermal runaway occurs at about 65 min.) 
During the long heat absorption period (0 min< t< 50 min), cap- 
turing the distributions in the three-dimensional model results in 
the lower rate of heat absorption from ambient at the surface and 
lower reaction rates at the inner core region. So, the thermal run- 
away appears to be delayed by a few minutes compared with the 
lumped model case. Due to the definite nonuniform distribution 
of temperature (the much higher than average temperature inside 
a cell), negative-solvent reaction progresses further on average 
so that the total heat release is greater to increase the maximum 
average cell temperature in the three-dimensional simulation 
case than in the lumped system analysis case. 

A sequence of component heat generation contours for the 
D50H90 size cell oven test simulation is presented in Fig. 12. The 
SEI decomposition reaction occurs much earlier (15—40 min) 
than the thermal runaway. Therefore, only the heat releases 
of positive-solvent and negative-solvent reactions are presented 
here. 


3.2.2. Localized heating simulation 

The three-dimensional model was used to simulate a local- 
ized heating event (such as an internal short circuit) to investigate 
the thermal abuse reaction propagation inside a D50H90 cell. A 
three-dimensional model is essential for simulating a localized 
heating event that is intrinsically multidimensional. A half-cell 
model with a symmetry plane was introduced. Internal hot spots 
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Fig. 12. Simulated sequence of component heat generation contours for a D50H90 cell in a 155 °C oven test using three-dimensional model. 
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were modeled by releasing localized energy in a short time in a 
very small region inside the cell core. Boundary conditions were 
also localized. Each boundary section can have various bound- 
ary conditions independently. In the following case, a specified 
amount of energy (here equivalent to 15% of stored energy in 
the cell) is released in a small portion of core volume (0.5% 
of total jellyroll volume) for 1 s. Initial temperature of the cell 
was 55°C. The cell was assumed to be cooled at the side sur- 
face with a heat transfer coefficient h =5 W m~? K and ambient 
temperature Tamb = 35 °C. 

Fig. 13 shows the spatial distributions of temperature and heat 
generation, and the shape of a reaction surface at 6s after the 
initial local heat release at the middle of the cell winding. The 
temperature contour presented in Fig. 13(a) shows that the initial 
hot temperature region mainly expands in the azimuthal and 
longitudinal directions. The initial reaction front propagating 
directions are more clearly seen in heat generation contours in 
Fig. 13(b). 

The pattern of thermal abuse reaction propagation following 
the initial local spot heating at the middle of the cell core wind- 
ing was investigated. The temporal variations of cell average 
temperature and total heat release are plotted in Fig. 14. Dur- 
ing the propagation of the abuse reactions through the cell, we 
observed several propagating stages where the cell average heat 
release rate increases suddenly. The corresponding temperature 
contours were produced and shown in Fig. 14. 
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Fig. 13. Computed spatial distributions of temperature and heat generation at 
6s after the initial local heat release. 
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(1) Initial propagation: In the beginning, the thermal abuse 
reactions propagate in the azimuthal and longitudinal direc- 
tions from the reaction ignition spot so that the reaction zone 
forms a hollow cylindrical shape. 

(2) Reaction at the core axis: Since the internal temperature dis- 
tribution largely depends on conduction heat transfer at this 
stage, the highest temperature appears at geometric center 
of the cell. As reactions and heat transport go on, the center 
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Fig. 14. Simulated pattern of reaction propagation due to local heating. 
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Fig. 15. Simulated iso-surface contours of SEI decomposition reaction propagation due to local heating. The shown numbers are time after the initial localized heat 


release. Each frame is 4s apart. 


of the cell starts to react and the reactions propagate along 
the cell axis where the nonreacted materials remain. 

(3) Reaction through the outer rim: As the cell is continuously 
heated, the outer rim of the cell is taken to the reaction from 
bottom to top. Since the container can is thermally connected 
to the bottom of the jellyroll, the bottom section of the cell 
is close to the thermal center, and heats sooner than the top 
section of the cell. 


The thermal abuse reaction propagation is more clearly seen 
at the consecutive iso-surface contours by presenting a sequence 
of SEI decomposition reaction completion surface (see Fig. 15). 


4. Conclusion 


Chemical reactions at elevated temperatures in Li-ion 
batteries were listed and formulated to be included in a three- 
dimensional Li-ion battery thermal abuse model. Reaction 
chemistry was implemented into a three-dimensional, finite- 
volume method battery model that addresses the effects of 
nonuniform distributions, thermal and electrical path design 
inside cells, localized heating and cooling, and cell component 
geometry and materials. Oven simulation indicates that cell size 
(or more precisely, heat transfer area per unit volume) greatly 
affects the thermal behavior of a cell in an oven test. Nonuni- 
form distribution effects were captured in the three-dimensional 
battery simulation. A localized heat release (representing a 
short circuit) was simulated, and the results were visualized 
to show the cell-internal propagation of the abuse reactions. 
The reactions initially propagate in the azimuthal and longitudi- 
nal directions to form a hollow cylinder-shaped reaction zone. 
Then the heat accumulation causes the center core of the cell 
to react. Finally, the reaction propagates to the outer rim of the 
cell. 

This study describes the development of a three-dimensional- 
capable thermal abuse model for lithium-ion cells and success- 


fully demonstrates multidimensional behaviors of thermally 
abused cells. Future work will include expanding the model 
to address various chemistries and cell-to-cell abuse reaction 
propagation in a multistring module. 
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